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Abstract 

We present a first-principles study on the interface between perovskite ferroelectrics (PbTiOa) 
and conducting magnetic manganites (Lai_ x Sr x .Mn03). We show that by switching the ferroelec- 
tric polarization, additional carriers are accumulated or depleted at the interfacial region of the 
manganite and that this change in carrier density can modify the magnetic spin configuration of 
the interfacial Mn, which is consistent with the experimentally observed anomalously large change 
in the magnetization. We also describe an unexpected purely interfacial phenomenon whereby the 
ferroelectric polarization of the interfacial region changes the magnetic energetics — a degree of 
freedom not present in bulk manganites. Theoretically, we show the the ground-state magnetic 
structure depends sensitively on the precise choice of Hubbard U parameter within the widely- 
used DFT+[7 class of exchange correlation functionals. We provide a simple Ising-like model that 
explains the evolution of the magnetic structure with U in tandem with a discussion of various 
different ways in which one might try to choose an appropriate U parameter. 
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I. INTRODUCTION 



Multiferroics have been one of the most intensively studied materials during the past 
decade jl, 2]. The coexistence of more than one order parameter in a single phase and their 
coupling may open new routes to the next generation of electronic devices. For instance, 
the possibility of controlling magnetization via external electric fields may find promising 
applications in spintronics. The origin of magnetoelectric multiferroicity lies in a nonzero 
magnetoelectric coupling which may occur due to many different mechanisms (for recent 
reviews, see {3, 4]). The magnetic properties of an intrinsic bulk magnetoelectric, of which 
Cr 2 3 is a prototype, can be modulated by an external field through the change of the 
magnetic cations' displacement relative to anions [5]. Extrinsic magnetoelectric couplings 
are typically mediated by strain: in composites of piezomagnetic materials combined with 
electrostrictive materials, external fields modulate the electric polarization, as well as the 
shape of the piezoelectric. This change in turn induces strain of the magnetic components 
and modifies the magnetization in the magnetostrictive material [6J. 

However, although these bulk mechanisms are well understood, the magnitude of mag- 
netoelectric couplings in bulk materials is generally small [7] , impeding their applications in 
electronic devices. Moving away from bulk materials, artificial heterostructures such as in- 
terfaces are promising candidates for realizing or even engineering magnetoelectric couplings. 
Recently, a variety of mechanisms were proposed to induce magnetoelectric coupling at inter- 
faces. At ferromagnet /ferroelectric interfaces, the interfacial bond length can be altered by 
ihe presence of ferroelectric polarization, for example in Fe/BaTiC>3 [8], Co2MnSi/BaTi03 
9] and FesC^/BaTiOs flp| . Another mechanism is to apply an external field and induce 
magnetization mediated by free screening carriers accumulated at ferromagnetic/dielectric 
interfaces, for example SrRuOs/SrTiOs ljj. A more complex but intriguing mechanism 



was recently described theoretically and experimentally [12J, [l3| at ferroelectric/conducting 
magnetic manganite interfaces. Our work focuses on this last class of systems, because the 
coupling between ferroelectric polarization and magnetization not only is of great impor- 
tance in fundemental sciences, but also finds very promising applications in memory devices 

At a ferroelectric/manganite interface, the presence of the ferroelectric polarization causes 
screening charges to appear at the interface due to accumulation or depletion of carriers in the 
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interfacial region. Much like SrRuOs/SrTiOs interfaces, the magnetization of the interfacial 
atoms can be enhanced due to the modification of carrier density around the interface 
because the magnetic moment of the atoms depends on the doping level (provided that the 
manganite is in the ferromagnetic phase). What is different in the ferroelectric/manganite 
system is that the accumulation of carriers not only changes the magnetic moment but 
can lead to an interfacial ferromagnetic-to-antiferromagnetic transition which reverses the 
directions of the moments and thus to a much larger magnetoelectric coupling. The work 
of Ref. [12I studied a representative heterostructure: BaTiOs/Lai-xBa^MnOs with x — 0.5. 
For bulk Lai-^Sr^MnOa, x = 0.5 is at the critical doping level separating ferromagnetic and 
antiferromagnetic phases, so the system is highly susceptible to magnetic changes with small 
changes of doping. What was found is that when the ferroelectric polarization is flipped, the 
magnetic moment of the Mn atoms in the second unit cell away from the interface is reversed. 
We note that the result is somewhat unintuitive as the carrier doping density is highest 
in the first layer at the interface which is most susceptible to change of magnetic phase. 
Experiments on the Pb(Ti .8Zr .2)O3/Lai_ ;c Sr :c MnO3 are performed for x = 0.2 [3] which 
is quite far from the boundary: a priori it is not clear whether enough screening charges can 
accumulate to drive the system over the magnetic phase transition. The experiments find a 
large magneto-electric coupling which is interpreted to originate from a spin-flip in the first 
unit cell of the manganite closest to the interface [131 ] . 

In this work, we comprehensively study this proposed magnetoelectric phenomenon 
at the ferroelectric/manganite interface using first-principles calculations. We choose 
PbTiOs/Lai-zSr^MnOs as our prototype for computational convenience. Although the ex- 
periments use Pb(Ti .8Zr 2 )0 3 as the ferroelectric [13| in order to reduce leakage currents, 
the key role of the ferroelectric is to induce screening charges due to its surface polarization, 
and we believe the compositional details are secondary. For the Lai_ x Sr x Mn03 thin film, 
we choose the nominal hole doping to be x = 0.2 (unless otherwise specified) in order to 
directly compare to experiment 131 ] . 

The original contributions of this paper are the following, (i) We are able to show that 
even though the manganite film has a low nominal doping of x = 0.2, the ferroelectric field 
effects can sufficiently change the carrier density to induce a magnetic transition at the inter- 
face. However, the precise ground-state magnetic structure of the interface (e.g., no spin-flip, 
spin-flip in first Mn layer, spin-flip in second Mn layer) depends sensitively on the choice 
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of exchange correlation functional and specifically the value of the Hubbard U parameter 
in the DFT+U approach, (ii) We clarify the physical reasons of the sensitivity based on a 
simple Ising-like nearest-neighbor model using bulk-derived parameters that well describes 
the computed dependences, (iii) We describe an unexpected purely interfacial effect that 
significantly modifies the energies of magnetic states at the interface: the ferroelectric polar- 
ization propagates into the first few manganite layers and the resulting rumplings of atomic 
layers are responsible for the modifications. This degree of freedom is not present in the 
bulk and potentially represents a further degree of freedom that can be exploited to modify 
and engineer material properties at oxide interfaces, (iv) We show that various approaches 
to find an appropriate U produce significantly different U values, some of which do poorly 
when used to compute results that compare to bulk properties of manganites; empiricism 
in the choice of U is an unfortunate necessity for manganites using state-of-the-art ab initio 
methods, (v) In the process of this work, we develop an alternative method for counting 
electrons on the Mn atoms in Lai-^Sr^MnOa thin films that is directly based on the electron 
density instead of the standard and widely used method based on projecting onto Lowdin 



or atomic-like orbitals [l6j. With this method, we can quantify the carrier distribution as a 
function of ferroelectric polarization and calculate the layer-resolved effective hole doping. 
The method is generally applicable to half-metallic oxide films. 

The remainder of the paper is organized as follows. We discuss computational details 
in Section [III We first study the magnetic phase transition of bulk Lai-^Sr^MnOs in Sec- 
tion [TTlAl A discussion of Hubbard U for bulk Lai_ x Sr x Mn03 is presented in Section UlI Bl 
The charge modulation at the interface is studied in Section IIVBI and Section IIVCI is de- 
voted to the discussion of spin modulation at the interface. We conclude in Section [V] A 
number of appendices contain further technical details. 

II. COMPUTATIONAL DETAILS 

Our calculations are performed using density functional theory within the ab initio su- 



jercell plane-wave approach 



171 ]. with the code PWscf in the Quantum-ESPRESSO package 



We employ ultrasoft pseudopotentials 18[. The semicore states and reference config- 



uration of each element are shown in Table [H We use the local spin density approximation 
(LSDA) [19] for the exchange correlation functional as well as the Hubbard U correction 
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TABLE I: The semicore states and reference configurations of our pseudopotentials. The cut-off 
radii are in units of Bohr. 



Atom 


Reference valence states 


< 


r p c 


r d 
' c 


Pb 


5d 10 6s 2 6p 2 


2.5 


2.5 


2.3 


Ti 


o 9o (in ll i 9 

3s ji 3p°3d 1 4s' i 


1.8 


1.8 


1.8 


Sr 


4s 2 4p 6 5s 2 


2.0 


1.8 


2.0 


La 


5s 2 5p 6 5d 1 6s 1 - 5 6p 05 


2.2 


2.0 


2.2 


Mn 


3s 2 3p 6 3d 5 4s 2 


2.0 


2.0 


2.0 


Pt 


5d 9 6s 1 6p° 


1.0 


1.0 


1.2 





2s 2 2p* 


1.3 


1.3 





method (LSDA+C) 20j to account for some of the strong electronic correlations on the 
localized d orbitals of Mn atoms. The plane wave basis energy cutoff and charge cutoff are 
35 Ry and 280 Ry, respectively. We use a Gaussian smearing width of 5 mRy when sampling 
the Brillouin zone. For bulk Lai_ x Sr x Mn03, the fc-grid sampling of the Brillouin zone is 
20 x 20 x 20 per formula unit. For interface calculations, the fc-grid sampling is 20 x 20 x 2 
where the z-axis is orthogonal to the interface. For variable cell relaxations, the convergence 
threshold for pressure is 0.5 Kbar. For atom relaxations, the convergence threshold for ev- 
ery force component is 26 meV/A. We have checked the convergence in total energies and 
structural parameters by further increasing the fc-point sampling and reducing the stress 
and force threshold, and observe no significant differences in key physical observables. 
The A-site Lai_ x Sr x alloying is treated by the virtual crystal approximation 



21 



22]. 



Appendix |X] describes tests on the accuracy of the virtual crystal approximation for our 
system: the results are highly satisfactory and consistent with earlier observations 23j 
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FIG. 1: (Color online) a) 20-atom unit cell of Pnma Lai-^Sr^MnOa. b) Experimental phase dia- 
gram of Lai-^Sr^MnOa. 'CI' means canted insulating. 'FI' means ferromagnetic insulating. 'FM' 
means ferromagnetic metallic. 'AFM' means antiferromagnetic metallic. 'PI' means paramagnetic 
insulating. 'PM' means paramagnetic metallic. The phase boundary of the ferromagnetic to A- 
type antiferromagnetic phase transition at low temperatures is highlighted by the bold dashed line. 
The panel b) is reproduced with permission from Ref. [3]. Copyright 2006 Institute of Physics 
Publishing. 

III. BULK MANGANITES 

A. Phase transition of magnetic ordering 



Since, as we will show, the magnetic properties of PbTiC , 3/Lai_. 2: Sr. 2 .MnC)3 interface can 
be understood qualitatively in terms of bulk Lai_ x Sr x MnC>3, we start the discussion with the 
phase diagram of bulk Lai-^Sr^MnOa under different conditions. The parent compound of 



Lax_2.Srj.MnO3 is LaMn03 which is an _l-type antiferromagnetic Mott insulator 



24|. Bulk 



LaMnOa has strong Jahn- Teller and GdFeOs distortions with Pnma symmetry 25j and its 
primitive cell is of size c(2 x 2) x 2 in units of the cubic perovskite. The smallest unit cell of 
LaMnO"3 has four formula units (20 atoms), as is illustrated in Fig. [jji. Chemically doping 
LaMnO"3 with Sr induces holes on the Mn <i-orbitals, leading to conduction and various 




Hole doping x Hole doping x 



FIG. 2: (Color online) Phase diagram of bulk Pnma La1_.cSra.MnO3. AE is denned in Eq. ([T]). 
Above zero, the ground state is ferromagnetic and below zero it is _l-type antiferromagnetic. a) 
Hole doping dependence within LSDA. b) Hubbard U dependence, c) Strain dependence, d) 
Structural distortion dependence. 

magnetic orderings. Fig. [TJd shows the experimental phase diagram of bulk Lai_~Sr x Mn03. 
A ferromagnetic to A-type antiferromagnetic phase transition occurs around x =0.5 doping, 
which is highlighted by the bold dashed line in Fig. [Tb. For a random alloy distribution, 
we assume that bulk Lai_ x Sr x MnC>3 has the same symmetry as LaMnC>3 (Pnma). In DFT 
simulations, we replace La with the fictitious atom Lai_-Sr- in the virtual crystal approxi- 
mation and calculate the energy difference between ferromagnetic ordering (F) and A-type 
antiferromagnetic ordering (A) as a function of doping x. 

AE = E(A) — E(F) (1) 

In addition to the doping dependence x, we also study the effect of strain, structural distor- 
tions and Hubbard U on the magnetic transition of bulk Lai_~Sr~Mn03. 
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1. Doping dependence 



Fig. |2k shows a representative bulk phase diagram of Lai_ x Sr x Mn03 as a function of 
doping x. The calculation is performed in the LSDA (i.e. Hubbard U—0). DFT-LSDA 
reproduces the ferromagnetic-to-antiferromagnetic phase transition that is observed in ex- 
periment. The doping dependence can be understood as follows: the ferromagnetic ordering 
is stablized by the double exchange mechanism which relies on the hopping of the lone e g 
electrons among the Mn 3+ ions 26(. With the increasing hole doping x, the itinerant elec- 
trons (i.e. Mn 3+ ions) are drained and the hopping processes are suppressed. Therefore 
the double exchange mechanism becomes less operative and the ferromagnetic ordering gets 
more unstable as x is increased. We need to point out that in the experimental phase dia- 
gram for x < 0.1, Lai-^Sr^MnOs is a spin-canted insulator (a magnetic structure close to 
A-type ant if erroma gnet ic ordering but the magnetic moment is not exactly cancelled due to 
weak spin-canting) [27]. As x increases to 0.2, there is an insulating-to-metallic transition 
and the appearance of the ferromagnetic ordering which is not reproduced in the DFT-LSDA 
calculation, since the LSDA ground state is metallic in the whole doping range. Turning on 
the Hubbard U does not change the metallicity of ferromagnetic Lax_j.Sra.MnO3. We argue 
that even though LSDA or LSDA+t/ is not accurate enough to reproduce the spin-canted 
ground state at very low doping nor the insulating-to-metallic transition around x ~ 0.2, 
it correctly produces the metallic ferromagnetic to metallic antiferromagnetic phase transi- 
tion at larger x , which is the key to understanding the spin-modified magnetic structure of 
manganites in the presence of ferroelectrics. 



2. Hubbard U dependence 

Fig. [2b shows a comparison of LSDA and LSDA+f/ calculations for bulk Lai__SrzMn(_3. 
The ferromagnetic-to-antiferromagnetic phase transition is well reproduced in both LSDA 
and LSDA+t/ calculations, but the transition point, i.e. the critical hole density where the 
ground state changes magnetic ordering, depends on the value of the Hubbard U. With an 
increasing U, the transition point moves to larger doping values while the overall features 
of the transition remain unchanged. The Hubbard U dependence originates as follows: 
antiferromagnetism is favored by the superexchange mechanism which involves the virtual 
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hopping of electrons between low and high energy sites with the same spin |28||. A larger 
U increases the energy splitting and thus the virtual hopping is suppressed. Therefore the 
superexchange mechanism is suppressed as U increases, and the antiferromagnetic ordering 
accordingly becomes less stable, resulting in the upward shift of phase transition curve 
(favoring ferromagnetism) . Empirically in order to correctly locate the transition point at 
the experimental value of x ~ 0.5, we need a Hubbard U in the range of 1 eV< U <2 eV in 
the LSDA+[7 approximation (as illustrated in Fig. [2b). 

3. Strain dependence 

Since the Lai_ x Sr x Mn03 thin film is grown coherently on an SrTiOs substrate, we also 
study the phase diagram of SrTiOs-strained Lai-xSr^MnOa and compare it with unstrained 
bulk Lai__ :r Sr :r Mn03 in Fig. [2b. Lai-zSr^MnOa in the whole doping range is under weak 
tensile strain (within 1%) when on an SrTiOs substrate. Tensile (compressive) strain re- 
moves the degeneracy of Mn e g orbitals and favors d x 2_ y 2 (d 3z 2_ T 2) orbitals due to the change 



of crystal field [29J. Based on the double exchange mechanism, ferromagnetism is isotropic 
with equal hoppings between Mn atoms along x, y and z directions. A-type antiferromag- 
netism is ferromagnetic in-plane and alternates its spin orientation layer by layer along the 



out-of-plane axis 30j. Due to the tensile strain, the occupancy of d 3d 2_ r 2 is lowered and 
the hopping between Mn atoms becomes essentially two-dimensional, suppressing ferromag- 
netism. Therefore with tensile strain, ferromagnetism is destablized and the whole transition 
curve is shifted downwards (favoring A-type antiferromagnetic ordering) as seen in Fig. [2b. 



4- Structural distortions 

Distortions away from cubic symmetry play a crucial role in the magnetism of manganites 



31] . Bulk Lai-zSr^MnOa has complicated structural distortions with Pnma symmetry (the 
unit cell is c(2 x 2) x 2 with 20 atoms). However, we also theoretically study 'artificial' 
Lax_j.Srj.MnO3 with only tetragonal distortions (the symmetry is P4 and the unit cell is 
1x1x2). The main reason we consider the high symmetry phase (-P4) and compare it to 
the low symmetry structure (Pnma) is computational: the P4 symmetry allows for the use 
of a smaller lxl interface unit cell which allows for simulation of much thicker films and 
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substrates. Therefore we need to understand the main differences, if any, between the two 
phases for what follows below. In addition, a comparison allows us to elucidate the role of 
structural distortions. 

Fig. [2U shows a representative phase diagram versus doping for both Pnma 
Lai_3.Sr3.MnO3 and P4 Lai_-Sr-Mn03 in the LSDA approximation. To understand these 
results, we begin with the fact that the effective hopping matrix element t between neigh- 



bouring Mn atoms depends on the Mn-O-Mn bond angle [32]. In the P4 case, the bond 
angle is 180° and the hopping is maximized, while in the Pnma case, the bond angle is 
smaller than 180° and the hopping is reduced. The double exchange mechanism depends 



ow- 



32]. 



linearly on this effective hopping matrix element t, while the superexchange mechanism 
ers the energy of antiferromagnetism by oc t 2 from second-order perturbation theory 
Therefore as we increase the hopping matrix element t, superexhange is more significantly 
enhanced than double exchange, thus favoring antiferromagnetism. Compared to the Pnma 
case, the P4 case has a larger effective hopping and the transition curve is shifted to favor 
antiferromagnetic ordering. This trend holds for both LSDA and LSDA+[7 calculations. 
Therefore, phenomenologically we can map Pnma La1_3.Sr3.MnO3 to P4 Lai_ x Sr x Mn03 by 
choosing an appropriate Hubbard U. We find that in order to reproduce the ferromagnetic- 
to-antiferromagnetic transition around x ~ 0.5 in the P4 case, we need 3 eV < U < 4 eV. 
By comparison, to locate the correct transition point for Pnma Lai_3,Sr3,Mn03, U must be 
in the range of 1 eV < U < 2 eV (see Fig. |2b). 



B. Choosing Hubbard U 



The DFT+U approach is commonly used to study manganites [25] . However, neither the 
choice of Hubbard U value nor the method of choosing it is unanimous. Obviously, one can 
choose U based on purely empirical considerations that use experimental data: for example, 
we showed above that when 1 eV < U < 2 eV, LSDA+[7 can correctly locate the experimen- 
tal critical doping density (x = 0.5) separating ferromagnetic and antiferromagnetic phases 
for bulk Pnma manganites. Below, we discuss two other reasonable-seeming methods one 
might consider to determine U. The approaches yield very different values of U that tend not 
to overlap and do not do well in comparison to experiment. In our opinion, unfortunately 
there is no reliable way to determine U in a theoretical a priori manner. Our opinion is that 
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a single-particle approach such as DFT+U will generally run into difficulties in describing 
strongly correlated system such as manganites, so that empiricism in choosing parameters 
is a necessary fact of life. Since the magnetic properties depend sensitively on the value of 
U, in our mind a more fruitful approach is to study a wide range of U to understand the 
trends versus U and especially why the trends take the form that they do instead of trying 
to make specific predictions based on some particular choice of U. (The U dependence of 
bulk manganite and ferroelectric/manganite interfaces are discussed in Sections IIII A 21 and 
IIV C\ respectively.) 



1. Bulk LaMn0 3 

First we may ask what U value properly describes the parent material: bulk LaMnOs. 
This value then may be a reasonable guess for the doped m ang anites. Taking into account 



the structural distortions by using a c(2 x 2) x 2 unit cell [25( and by relaxing all degrees 
of freedom, we calculate the total energies of different magnetic orderings and find their 
energy sequence as a function of U. The result is shown in Table [Til In particular, we 
explicitly list AE, defined by Eq. flTJ in the table. We can see that within a wide range 
of U, the ground state is not the experimentally observed A-type antiferromagnet, nor is 
there any tendency that ferromagnetism could yield to antiferromagnetism in the large U 
limit. However, reproducing the insulating properties of the A-type antiferromagnetic phase 
requires U > 4 eV. 

We note that one can perform self-consistent calculations on bulk LaMnC>3 using the 
experimental lattice parameters and atomic coordinates. It is possible to stabilize an insu- 
lating A-type ground-state for U < 1 eV, as shown in Table IHH For a comprehensive study 
of bulk LaMnO^ studied with a variety of exchange correlation functionals and basis sets, 



please refer to [33j and references therein. Unfortunately, the reproduction for the correct 
ground state when using experimental structures is not of great value for our study: we have 
a non bulk-like interfacial system where the in-plane lattice constants are fixed via epitaxy 
to a substrate and all remaining degrees of freedom must be relaxed, so we must return 
to Table [ITJ It would seem the best choice is either U = (which stabilizes the incorrect 
ground-state by the least energy) or U > 4 (which makes the v4-type phase insulating). As 
shown above, neither choice is satisfactory in reproducing the experimental x ~ 0.5 phase 
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boundary for the doped manganites. 



TABLE II: LSDA+C/ study of bulk LaMn0 3 . The italic M means metallic and the itallic / 



means insulating. F refers to ferromagnetic ordering. A, C and G refer to ^4-type, C-type and 

□ 

G-type antiferromagnetic ordering, respectively [30]. AE is the energy difference between the 
ferromagnetic ordering and A- type antiferromagentic ordering per Mn atom, defined by Eq. (T1J. 
The unit cell is orthorombic. The experimental value of lattice constants is: a = 5.742 A, b = 7.668 
A and c = 5.532 A 



251 ] . The calculated lattice constants are for A- type antiferromagnetic ordering 
because the experimental ground state is ^4-type antiferromagnetic. In the parenthesis lists the 
relative difference between experimental and theoretical lattice constants. 



U (eV) 


magnetic ordering 


AE (meV) 


a(A) 


6(A) 


c(A) 





F(M) < A(M) < C(M) < G{M) 


40 


5.402 (-5.9%) 


7.468 (-2.6%) 


5.458 (-1.3%) 


2 


F(M) < A(M) < C{M) < G(I) 


61 


5.567 (-3.0%) 


7.560 (-1.4%) 


5.435 (-1.8%) 


4 


F(M) < A(I) < C(I) < G(I) 


65 


5.644 (-1.7%) 


7.584 (-1.1%) 


5.448 (-1.5%) 


6 


F(M) < A(I) < C(I) < G(I) 


76 


5.699 (-0.8%) 


7.624 (-0.6%) 


5.465 (-1.2%) 


8 


F(M) < A(I) < C(I) < G(I) 


98 


5.743 (0.02%) 


7.694 (0.3%) 


5.482 (-0.9%) 



2. Linear response approach of self- consistent U 

Second, we may ask for a purely ab initio approach that delivers a value of U appropriate 
for the s yste m within the framework of DFT+U itself. This is the linear response approach 



of Refs. 



35|. We focus on P4 Lai^Sr^MnOs (x = 0.2) as an example. We run a series 



of linear response calculations 34] on 2 x 2 x 2 unit cells of P4 Lai_ x Sr x Mn03. The ground 
state is calculated using LSDA+f/ with a range of < U in < 5 eV. For each value of U in , 

n 

we use the extrapolation scheme in Ref. [34| to get the converged value of U out . Then we 
collect all the converged U out as a function of Ui n and extract out U sc f [35] from the linear 
region. Our final value is U sc f = 5.8 eV. As discussed above, to reproduce the experimental 
x ~ 0.5 boundary for PA Lai_ x Sr x Mn03, we require 3 eV< U < 4 eV. The self-consistent 
U is significantly higher. 
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TABLE III: LSDA+C7 study of bulk LaMn03 using experimental coordinates and lattice constants. 
The italic M means metallic and the itallic / means insulating. AE is the energy difference between 
the ferromagnetic ordering and A- type antiferromagentic ordering per Mn atom, defined by Eq. ([!]). 



U (eV) 


magnetic ordering 


AE (meV) 





A(I) < F(M) 


-15 


1 


A(I) < F(M) 


-4 


2 


F{M) < A(I) 


4 


3 


F{M) < A(I) 


10 


4 


F{M) < A(I) 
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IV. FERROELECTRIC /MANGANITE INTERFACES 
A. Methodology 

Our computational supercell for interface calculations is schematically illustrated in 
Fig. [3j The x and y directions of the simulation cell are subject to periodic bound- 
ary conditions and their lengths are fixed to our computed theoretical lattice constant of 
SrTiC>3 a = 3.85 A (1.5% smaller than the experimental value), because in experiments 



Lai_ x Sr x Mn03 is epitaxially grown on a SrTiOs substrate 13J]. In order to directly com- 
pare to the experiment, all the results shown below are from calculations with the nominal 
doping level chosen as x = 0.2 (unless otherwise specified). In addition to PbTiC>3 and 
Lai_3.Sr3.MnO3, we also include the electrode Pt to provide an electron reservoir and ~ 
20 A vacuum to separate periodic copies of the slabs. We strain the in-plane lattice constant 
of the entire slab structure to that of bulk SrTiOs to impose the epitaxial strain from the 
substrate. In the simulation cell (Fig. [3]) and in most of our calculations, we do not include 
a SrTiOs substrate explicitly in order to keep the computations from becoming unwieldy in 
scale. However, in Appendix |B] we present a few calculations that do include the SrTiOs 
substrate explicitly, and it is shown that the interfacial structural and magnetic properties 
between PbTiOs/Lai^S^MnOs are well converged when the Lai_3;Sr3;Mn03 film is 4 unit 
cells or thicker. In addition to reducing the computational burden, the absence of a SrTiOs 
substrate creates a manganite surface that allows us to apply a hole counting method which 

13 



PbO|Mn0 2 



a 



SrTiCK 




Pt 



Pb 



Ti Mn La/Sr O 



FIG. 3: (Color online) Illustration of the computational supercell. The dashed line highlights the 
PbO/Mn02 interface. The whole structure is coherently strained to the lattice constant of SrTiO"3. 
Vacuum ~ 20 A thick is introduced to separate periodic copies. An XO (X=La/Sr) atomic layer 
faces the vacuum. 

can much more accurately calculate the hole spatial distribution than the use of Lowdin 
orbitals |16| (see Appendix ICl for details). Since the SrliC^ substrate is typically TiCV 
terminated and the manganites are in principle deposited stoichiometrically and epitaxially 



on the SrTiOs, the resulting ferroe 
Using the Berry phase method 



lectric/manganite interface is taken to be PbO/Mn02- 
3^], we find that SrTiCvstrained PbTiOs has bulk po- 
larization 0.74 C/m 2 . For the two different directions of ferroelectric polarization, we define 
two distinct states: the accumulation state in which extra holes are induced into the in- 
terfacial Lai^Sr^MnOs and depletion state in which extra electrons are induced into the 
interfacial Lai-zSr^MnOs (i.e., holes are driven out). One unit cell of PbTiOs in the interior 
is fixed to the bulk ferroelectric PbTiOa positions, an choice that simulates the behavior of 
a thick PbTiOs film. All remaining atomic coordinates in the slab are relaxed. We need to 
mention that our choice of boundary condition on ferroelectrics is consistent with the exper- 
iment 13] in which a thick film of 250 nm Pb(Zro.2Tio.s)03 is deposited on Lai_ x Sr x Mn03. 
However, the boundary condition on the ferroelectric could be different, depending on the 
experiments to be studied. For example, in Ref. [15], three unit cells of BaTi03 adjacent 
to a Lao.7Sr .3Mn03 film are fully relaxed without the presence of Pt electron reservior, in 
order to simulate the ultra-thin ferroelectrics used in other experiments |l4|. Interestingly, 
the results of magnetoelectric coupling in Ref. [15j] are consistent with ours, described below. 
In Fig. HI we show the cation-oxygen z-axis displacements of a representative PbTiC>3 
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FIG. 4: (Color online) Cation-oxygen displacements along z-axis for the PbTiO"3 film inside the 
Pt/PbTi03/Lai_ x Sr x Mn03 heterostructures. The two purple solid lines highlights the one unit 
cell of PbTiO"3 that is fixed to bulk positions. The two orange dashed lines show the interfaces: 
left is the one facing Pt electrodes and right is the one that faces Lai_ x Sr x .Mn03. 



thin film within the slab structure. The single fixed unit cell is highlighted by the two 
solid purple lines. The two interfaces (one faces La 1 _ x Sr 2 .Mn03 and the other faces Pt 
electrodes) are shown by the orange dashed lines. The sign of the displacements indicates 
polarization directions. We can see that in the relaxed PbTiC>3, there is no reversal of 
ferroelectric polarization and the magnitude of polarization is homogeneous. In addition, 
we also calculate an artificial state in which one unit cell of PbTiC>3 is fixed to be paraelectric 
(i.e., zero cation-oxygen rumpling in the (100) atomic plane of the fixed unit cell). 

Finally, we mention that most of the results presented below are calculated for an in- 
plane c(2 x 2) unit cell which is compatible with the structural distortions found in bulk 
Pnma Lai_ x Sr x Mn03. Such calculations are referred to as c(2 x 2) Lai_ x Sr x Mn03 inter- 
face calculations. In order to converge the hole distribution versus Lai-^Sr^MnOs thickness 
without inordinately increasing the computational burden, we increase the thickness of man- 
ganites by reducing the in-plane cell to 1 x 1. Those calculations are referred to as 1 x 1 
Lai_ x Sr x Mn03 interface calculations. 
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FIG. 5: (Color online) Left panel: hole distributions of the accumulation and depletion states of 
c(2 x 2) Lai_ x Sr x Mn03 (U = 1 eV). Middle panel: Hubbard U dependence of hole distribution in 
the accumulation state. Right panel: Hubbard U dependence of hole distribution in the depletion 
state. 

B. Charge modulation 

We first study the effect of charge modulation from switching the ferroelectric polarization 
of PbTi03. In the presence of ferroelectric PbTiOs, the charge density of Lai-^Sr^MnOa at 
the interface differs from its bulk value because the polarization of PbTi03 terminates at the 
interface and results in the surface charge (the surface charge density is a = Pn = P z ). Since 
Lax_j.Srj.MnO3 is metallic, this surface charge induces screening charge in the Lai_ x Sr x Mn03 
equal in magnitude but opposite in sign to the surface charge. When the PbTi03 switches 
its polarization, the surface charge changes sign and so does the screening charge. Therefore 
a net change of charge density (Act = 2P Z ) is induced in the Lai_ x Sr x Mn03 thin film. 

The left panel of Fig. [5] shows the hole distribution of a 4 unit cells thick 
Lax_j.Srj.MnO3 film on PbTi03. The method of counting holes is described in the Ap- 
pendix O The nominal doping is x = 0.2, highlighted by the dashed line. This calculation 
is performed on c(2 x 2) Lai_j,Sr x Mn0 3 with LSDA+f/ (U = 1 eV). The PbO/Mn0 2 in- 
terface is at layer 1, and layer 4 is the artificial surface. As expected, the hole distribution 
accumulates (depletes) at the interface when the surface charge is negative (positive). 

The middle and right panel of Fig. |5] show the Hubbard U dependence of spatial hole 
distribution for the accumulation and depletion states, respectively. The calculation is per- 
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formed on the same structure as in the left panel of Fig. [5j Since the induced holes (or 
electrons) extend into the Lai-^Sr^MnOs within the screening length (which does not 
strongly depend on the correlation), it is not surprising that Hubbard U does not signif- 
icantly change the hole distribution. As a good approximation, we assume that the hole 
distribution does not depend on Hubbard U. 

In order to get a spatial distribution of holes that is well converged with the manganite 
thickness, we run a calculation with 8 unit cells of 1 x 1 Lai^Sr^MnOs . This calculation 
is performed with LSDA+[7 (£7=4 eV). The results are shown in Fig. [6j From Fig. [6^, 
the screening length of the accumulation state is estimated to be 3 unit cells while the 
screening length of depletion state seems to be only 1 unit cell. The asymmetry is due 
to the fact that there are two factors affecting the hole distribution. One is the induced 
screening charge and the other is the presence of PbO /MnC>2 interface itself. To demonstrate 
the role of the interface, we perform a test calculation in which the PbTiC>3 thin film is 
forced to be paraelectric and find (see the orange triangle symbols in Fig. EK) that the 
resulting hole distribution is not uniform nor equal to the nominial doping (x = 0.2) at 
the interface. This non-uniform hole distribution can be considered as a background, owing 
to the chemistry of the PbO/Mn02 interface. If we substract the hole distributions of 
accumulation and depletion states from this background, we can see that the "net" hole 
distributions of accumulation and depletion states now become more symmetric, with the 
screening length of depletion state a little larger than that of accumulation state (Fig. [6b)- 
This is consistent with the Thomas-Fermi picture that depletion states have less carriers 
(holes) and therefore a larger screening length. To further verify this Thomas-Fermi picture, 
we perform the same calculation with a layer nominal hole doping x = 0.5 and find very 
similar results (see Fig. [6fc and d) . 

C. Magnetization modulation 

In this section, we study in detail whether the charge modulation can induce a spin 
modified configuration in the ground state. This means that the change of the magnetization 
is not simply proportional to that of the charge density (i.e. simple filling/emp tying of Mn 
orbitals with fixed spin polarization) but involves a more dramatic change of magnetic 
structure at the interface. The mechanism is as follows: in the accumulation state, the 
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FIG. 6: (Color online) a) Hole distributions of 1 x 1 Lai-^Sr^MnOs 8 unit cells thick with x = 0.2 
(U = 4 eV). In addition to accumulation and depletion states, a test calculation in which PbTiOa 
is fixed to be paraelectric is done and the resulting hole distribution is shown with orange triangles, 
b) "Net" hole distributions of 1 x 1 Lai^Sr^MnOs 8 unit cells thick with x = 0.2 (U = 4 eV): the 
paraelectric state background (see a)) is substracted from the hole distributions of accumulation 
and depletion states, c) Hole distributions of 1 x 1 Lai_ a Sr x Mn03 8 unit cells thick with x = 0.5 
(U = 4 eV). d) "Net" hole distributions in 1 x 1 Lai^Sr^MnOs 8 unit cells thick with x = 0.5 (U 
= 4 eV). 

local hole distribution adjacent to the interface could be higher than x = 0.5, the critical 
value for the ferro-to-antiferromagnetic transition. Therefore the spins at the interfacial 
region could flip. However, in the depletion state, such a local spin-flip is not expected to 
occur. Therefore from now on, we only focus on the accumulations state. In order to study 
whether this local phase transition does occur at the interface by switching the polarization, 
we consider three relevant spin configurations (F, Al and A2), illustrated in Fig [71 When 
all the spins are ferromagnetically coupled, this configuration is denoted as F (Fig. [7Ji). If 
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FIG. 7: (Color online) Illustration of different magnetic structures, a) Ferromagnetic configuration, 
denoted by F. b) Spin flips in the first unit cell, denoted by Al. c) Spin flips in the second unit 
cell, denoted by Al. 

the spin is nipped in the first unit cell of manganite from the interface, this configuration is 
denoted by Al (Fig. [7}d). Finally, if the spin is flipped in the second unit cell of manganite, 
then we denote it by Al (Fig. [7b). We address three important and related questions below: 
i) whether the ground-state magnetic structure depends on Ul ii) given a reasonable U, 
whether the manganite nominal doping x could change the final magnetic structure? iii) 
how the structural distortions at the interface may affect the magnetic structures? 

1. Hubbard U dependence 

We obtain the total energies of these three spin configurations with a range of Hubbard 
U and collect all the results in Fig. [HJ We use the following definitions of energy differences: 

A£i = E(A1) - E(F) (2) 

AE 2 = E(A2) - E(F) (3) 

From Fig. |8]we can see that the magnetic structure of the ground state evolves with Hubbard 
U. When U is small {U < 0.9 eV), the ground state has the magnetic structure of Al. With 
V increasing (0.9 eV < U < 2.4 eV), the ground state evolves into Al. When U gets even 
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FIG. 8: (Color online) Energy sequence dependence on Hubbard U. AEi = E(A1) — E(F) and 
AE2 = E(A2) — E(F). The labels (F, Al and A2) show the ground state in different regions of 
Hubbard U. The boundary between F and Al is 2.4 eV and that between Al and A2 is 0.9 eV. 

larger (U > 2.4 eV), we have F as the ground state. Therefore any prediction of the magnetic 
properties of the interface depends greatly on the choice of U. Before we pick a reasonable 
value of U, we need to understand why the magnetic structure is so sensitive to the Hubbard 
U. The answer is that U changes the bulk phase transition point so that for the same hole 
distribution, the preferred local magnetic phase also changes. We can see from Fig. [8] that 
increasing the Hubbard U drives the local phase at the interface from antiferromagnetic to 
ferromagnetic, which is consistent with the U dependence in the bulk (see Fig. [2)3). In order 
to more quantitatively describe the energy sequence, we construct an Ising-like model which 
is based on the interaction between nearest neighbor Mn magnetic moments 39(: 



E 



} J JijUiirrij 

<ij> 



(4) 



where < ij > range is over all nearest neighbors and Wj is the magnetization in each Mn02 
layer of the manganites. The labelling of manganite layers is shown in Fig. [71 We assume 
that the hole spatial distribution does not sensitively depend on magnetic structures [401 ]. 
and obtain: 
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(5) 



AE 2 = 2Ji 2 |mim 2 | + 2J 2 3\m 2 m 3 \ (6) 

In order to get an energy sequence, we need to know the signs of J\ 2 and J23. From the 
bulk calculations, at a given hole doping x and assuming half-metallicity, the magnetization 
is related to the hole doping x by m = (4 — x)/j,b where /i# is the Bohr magneton. The 
exchange coupling J can be extracted out by: 

_ E(A) - E(F) 
J ~ 2m 2 [0 

where the energy difference AE = E(A) — E(F) is from the bulk calculations, shown in 

Fig. El J changes sign at the transition point. From Eq. (EJ), J is positive for ferromagnetic 

phase and negative for A-type antiferromagnetic phase. At the interface, however, the hole 

spatial distribution is not uniform (see Fig. E]). We assume that the interface coupling 

is that of bulk La 1 _. r Sr. r Mn03 but for a doping value that is the average of the neighboring 

layers i and j: 

- J bu i k f^i] (8) 



2 

We need a final good approximation, which is verified in Fig. [51 that the hole distribution 
does not sensitively depend on Hubbard U. Based on Eq. ([ME]); we start with a large Hub- 
bard U. Since large U favors ferromagnetism (see Fig. [2b), the bulk phase is ferromagnetic 
and both J 12 and J23 are positive. Thus < AE\ < AE 2 and we have the following energy 
sequence: F < Al < A2. We denote this by case 1. With a decreasing U, the transition 
point is moved to smaller hole doping region. Noting that the hole distribution monoton- 
ically decays from the interface (see Fig. [6]), we always have (x\ + x 2 )/2 > (x 2 + x 3 )/2. 
Hence J\ 2 changes sign earlier than J23 as U decreases. If U is in such a range that J\ 2 just 
becomes negative but J23 > 0, we have AE\ < and AE 2 > 0. The energy sequence is now 
Al < F < A2, which is denoted by case 2. As U further decreases, so that J\ 2 becomes very 
negative and J23 remains positive but J^lwiil + </23| m 3| < 0, then we have AE\ < AE 2 < 0. 
The energy sequence becomes Al < A2 < F. This is case 3. With U further decreasing, the 
bulk phase becomes always antiferromagnetic, both J 12 and J 2 3 become negative, and we 

21 



TABLE IV: The energy sequence predicted from the simple model and the comparison with the 
DFT calculations with different Hubbard U. 



model 


DFT 


case 


energy sequence 


U (eV) 


AE X (meV) 


AE 2 (meV) 


1 


< AE\ < AE 2 


3 


14 


76 


2 


AE 1 < < AE 2 


2 


-8 


27 


3 


AEi < AE 2 < 


1 


-34 


-27 


4 


AE 2 < AEi < 





-70 


-120 



have AE 2 < AE% < 0. The final possible energy sequence is A2 < Al < F, which is denoted 
by case 4. These four energy sequences exhaust all the possibilities and are summarized in 
Table HVl Now we compare the DFT results (see Fig. [8]) to the energy sequence predicted 
from the model (see Table HV|) . As the Hubbard U evolves from to 4 eV, we find all four 
cases. For example, U = 3 eV corresponds to F < Al < A2\ U = 2 eV to Al < F < A2; 
U = 1 eV to Al < A2 < F and U = eV to A2 < Al < F. The exact boundaries of 
Hubbard U for each energy sequence can be found in Fig. [HJ 

Since the Hubbard U changes the transition point and the magnetic structure of the 
ground state of the PbTiQcj /Lai _. T Sr r MnO.s interface, we need to determine what is the 



reasonable value of U. Following Ref. 



12], we argue that because the magnetic struc- 



tures sensitively depend on the transition point, we need to choose a range of U so that 
the ferromagnetic-to-antiferromagnetic transition occurs around x ~ 0.5. From bulk cal- 
culations, we know that as 1 eV < U < 2 eV for Pnma Lai-^Sr^MnOs, this criterion is 
satisfied. On the other hand, when U is in this range, the magnetic structure of the ground 
state is always Al. Therefore by switching the PbTiC>3 polarization, we do find a spin- 
modified configuration in the DFT simulation, provided that our choice of U is reasonable. 
This prediction is consistent with the recent experiment [13j which observes an anomously 
large change in the magnetization as the polarization of ferroelectrics is switched and which 
assigns this to a spin-flip on the Mn atom closest to the interface. 
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FIG. 9: (Color online) Energy dependences of the interfacial magnetic states as a function of the 
doping x of the Lai-^Sr^MnOa. AE X = E(A1) - E(F) and AE 2 = E(A2) - E(F). The boundary 
between the phases is at x c ~ 0.22 denoted by the vertical dashed line. Al is stable to the left and 
A2 to the right of the boundary. The results are based on LSDA+C/ with U = 1 eV. 

2. Hole dependence 



The ground-state magnetic structure we found above (Al configuration) is consistent 
with the experimental conjecture, but it is different from the A2 configuration found using 
DFT+GGA for the similar multiferroelectric structure BaTiOs/La^Ba^MnOa with x = 
0.5 [12|| . We find that the reason for the differing ground state magnetic structure is due to 
the doping x dependence of the system. Specifically, for U = 1 eV, we calculate the energies 
of the Al and A2 interfacial states versus doping x and present the results in Fig. |9j 

When the nominal doping x is near the bulk magnetic transition point x = 0.5, the Fig. [9] 
shows that the ground state should be A2 which is consistent with Ref. 12]. However, 
the actual experimental doping x = 0.2 is far below x = 0.5, the ground state should be 
Al. We calculate the phase boundary between Al and A2 to be x c ~ 0.22 for U = 1 
eV. This transition from Al to A2 is easily understood in the framework of our Ising-like 
model. Each pair of neighboring Mn magnetic moments could be either ferromagnetically or 
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antiferromagnetically coupled depending on the number of holes on the pair of Mn atoms. 
When the doping x is low, the ferroelectric modulation of the hole density must drop to a 
low value within a few unit cells from the interface (see Fig. E]) and thus only the first two 
moments (m l5 m 2 ) will be paired antiferromagnetically, which leads to Al. However, when 
the doping x is high, the larger values of the hole density means that both the (mi, 777.2) an d 
(777.2,777.3) pairs couple antiferromagnetically, which leads to A2. As a final remark, we note 
that our Ising-like model predicts that the transition doping x c between Al and Al should 
depend on Hubbard U: since increasing U favors ferromagnetism, larger U will increase x c 
(i.e, require more holes for antiferromagnetism). 



3. Structural distortion dependence 

Due to the presence of ferroelectric polarization, significant distortions that deviate from 
bulk manganites are observed in the relaxed ground state structures. Concerning each 
oxygen octahedron that encloses Mn atoms, we calculate c/a ratio and rumplings S/a in 
each Mn02 layer, where c is the distance between the two apical oxygen atoms along the 
z direction, S is the vertical displacement between Mn and O, and a is the lattice constant 
of SrTi03 substrate. The results are summarized in Table |V] Since the spin-flipped process 
occurs at the interface in the accumulation state, we only show c/a ratio and S/a of the 
first and second manganite layers from the interface, and from now on the discussion is 
constrained to the accumulation state. From Table |V] we can see that in the accumulation 
state, there are sig nificant polar distortions at the interface (S/a is as large as 6%). It was 



shown in Ref. [12| that the spin-flipped process is mainly of electronic origin rather than 
due to the polar distortions at the interface. In this section, we use detailed comparisons 
to show that though the spin-modified configuration is due to electronic reconstructions, 
polar distortions need to be taken into account in order to make a quantitative (instead of 
qualitative) link between the interface phase and bulk phases. 

Now we look at the Ising-like model Eq. ([ME]) more closely. The model is based on the 
assumption that the local magnetic structure can be predicted from bulk manganites of the 
same hole doping. In Table Wl\ we list the energy difference between F and Al interfacial 
configurations from the interface calculations, defined by 
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TABLE V: The c/a ratio of each oxygen octahedron that encloses Mn atoms and rumplings 5/a 
of each MnC>2 layer for both accumulation and depletion states, c is the distance between the 
two apical oxygen atoms along the z direction. 5 is the rumplings of MnC>2 layer, and a is the 
lattice constant of SrTiC>3 substrate. Layerl and layer2 refer to the first and second unit cell of 
manganites from the interface. 



accumulation state 


depletion state 




layerl 


layer2 




layerl 


layer 2 


c/a 


0.97 


0.94 


c/a 


1.05 


0.98 


5/a 


0.06 


0.02 


5/a 


0.01 


0.001 



AEj = E(A1) - E(F) (9) 

where E(F) and E(A1) are the total energies of F and A\ configurations, respectively. 
We also calculate the average hole density between the first and second layers, i.e. x = 
(x\ + X2)/2. Next, we list the bulk energy difference AE B , denned as 

AE B = ^(E B (A)-E B (F)) (10) 

where E B (F) and E B (A) are the total energies of SrTiOs-strained Lai^Sr^MnOs with 
ferromagnetic and A-type antiferromagnetic ordering, respectively. The factor | is included 
because in the bulk form wherever a Mn atom flips its spin, there are two Mn-Mn bonds 
involved owing to periodic boundary conditions, whereas at the interface a Mn spin flip only 
involves one Mn-Mn bond. Hence, we need a factor | so that both AEj and AE B describe 
the energy difference per Mn-Mn bond. The nominal hole doping x is chosen as the same 
as x from the supercell calculations. Table |VT] shows that although the trend versus U is 
the same in both supercell and bulk calculations, the magnitudes of AE do not agree at all. 
There must be something at the interface which is absent in bulk phase and significantly 
affects the energy difference between antiferro- and ferromagnetism. 

We find that, due to the presence of ferroelectric PbTiC>3, strong polar distortions are 
induced at the interface layer of MnC>2 in the accumulation state (Table IV]) . as is illus- 
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TABLE VI: The comparison of energy difference between the interface calculations and bulk 
Lai_ x Sr x Mn03 calculations, x = (x\ + X2)/2 is the average hole in the first and second 
Lai_ x Sra;Mn03 layers closest to the interface. For different Hubbard U, x does not change signif- 
icantly. AEi is the energy difference between phases Al and F from the supercell calculations. 
AEb is the energy difference of SrTi03-strained bulk Lai-xSr^MnOs between A-type antiferro- 
magnetism (A) and ferromagnetism (F). AEbp is the energy difference between E(A) and E(F) 
of SrTi03-strained bulk Lai_ x Sr x .Mn03, with the interfacial polar distortions manually included 
and c-axis optimized. For bulk calculations (AEb and AEbp), the nominal hole density is chosen 
as x . 



U (eV) 


X 


AEi (meV) 


AE B (meV) 


AE B p (meV) 


1 


0.7 


-34 


-112 


-47 


2 


0.7 


-8 


-70 


-17 


3 


0.65 


14 


-24 


16 


4 


0.6 


42 


19 


57 



trated in Fig. ITUa . The cation-oxygen vertical displacement in MnC>2 at the interface is 
5 = 0.2 A. Such a strong polar distortion (distinguished from the structural distortions of 
Pnma symmetry) is absent in bulk Lai-^Sr^MnOs and is a direct consequence of the ferro- 
electric/manganite interface. In order to show explicitly how this interfacial polar distortion 
affects the energy difference between F and A-type magnetic orderings, we perform the fol- 
lowing thought experiment, which is schematically illustrated in Fig. [TUb . We focus on the 
two layers of Lai-^Sr^MnOa closest to the interface in interface calculations, use their relaxed 
atomic positions and choose an x that is equal to the average hole doping x — (x\ + x%)/2 
from interface calculations. In this way, we create such artificial La^^Sr^MnOa with the 
same built-in polar distortions and the same average hole doping as the two manganite layers 
at the interface. We impose periodic boundary conditions on this artificial Lai_ s S%Mn03, 
fix all atom positions and in-plane lattice constants, and optimize the c-axis to minimize 
the out-of-plane stress. We tune c separately for both ferromagnetism and A-type antifer- 
romagnetism. A similar energy difference AEbp (subscript P means 'polarized') is defined 
as 
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FIG. 10: (Color online) a) Strong polar distortion is induced at the interface layer of MnC-2 due to 
the presence of PbTiO"3. b) Schematics of copying the two interfacial Lai-^Sr^MnOs layers and 
forming artificial bulk Lai_jSrj-Mn03 with built-in polar distortions. 

AE BP = i (E BP (A) - E BP (F)) (11) 

where E BP (F) and E BP (A) are the total energies of the artificially constructed 
La!_^S%Mn03 with ferromagnetic and A-type antiferromagnetic ordering, respectively. We 
can see from Table [VI] that AE BP is much closer to AEj than the raw bulk data AE B , 
demonstrating that in order to quantitatively connect the phase evolution of the ferroelec- 
tric/manganite interface from the bulk manganite phases, the polar distortion induced in 
interfacial manganites is an essential ingredient in modelling. 



V. CONCLUSION 



We have presented a systemic study of the PbTiOa/Lai-^Sra-MnOs interface as a proto- 
type for ferroeletric/manganite interfaces. We are able to show that the screening charges 
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produced in the manganite in response to the ferroelectric surface charge are sufficient to 
change the magnetic state of the interfacial manganite from ferromagnetic to antiferromag- 
netic, in agreement with experimental observations and interpretations. In the process, we 
have developed a method to accurately count the layer-by- layer hole distribution in the 
manganite thin film which allows us to perform quantitative analysis of the system. For 
example, it allows us to create a simple Ising-like model of the interfacial magnetism that 
uses bulk parameters to reproduce the computed behaviors. 

One of main theoretical findings is that the ground-state magnetic state depends sen- 
sitively on the value of U chosen in the LSDA+[7 computation. We show that different 
reasonable-seeming approaches to determining U, and in particular some that are ab initio 
and deliver a U value appropriate to LSDA+[7 self-consistently, yield significantly different 
U values. Not all the values do well when compared to experiment. By asking that the 
LSDA+£7 calculation should correctly reproduce the critical hole doping density separating 
the ferromagnetic and antiferromagnetic phases, we are able to find a narrow range of U 
values that also produce a straightforward interfacial magnetic ground-state structure (the 
Al configuration) whereby the manganite layer with the highest doping has the strongest 
magnetic response. Clearly, our conclusions on the magnetic ground-state are not ab initio 
as they involve significant experimental input. In our opinion, the unsatisfactory situa- 
tion vis a vis choosing the U value is due to the limitations of the single-particle DFT+£7 
method itself when applied to a complex and strongly correlated electronic system such as 
manganites: the theory is not accurate enough for the material, so some level of empiricism 
is unfortunately necessary. 
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Appendix A: Tests of the virtual crystal approximation 

We treat the A-site Lai-^Sr^ alloying in Lai_ x Sr x Mn03 with the virtual crystal approx- 



imation 
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22]. This approximation involves replacing the two elements by a fictitious one 
whose electron number is (1 — x)Nj j3i + xNs t , where and iVg r are the number of electrons 
of the La and Sr pseudo atoms, respectively. We stress that i) since the magnetic properties 
originate from Mn <i-electrons and v4-site atoms serve to donate electrons, we expect that 
the virtual crystal approximation is reasonable to describe the magnetic phase transition 
of manganites in the random distribution; ii) since the chemical properties mainly depend 
on the valence electrons, the approximation we make here is expected to be also good for 
Ca and Ba, and iii) our choice of pseudo potential and valence electrons shown in Table |T] 
ensures a very smooth interpolation between La and Sr as their pseudo valence electrons and 
nuclear pseudo charges only differ by one elementary charge. We perform simple tests of 1:1 
Sr-La alloying (i.e. 50% alloying) in a c(2 x 2) x 2 unit cell and the results are compared 
to the x = 0.5 virtual crystal approximation as shown in Table IVlTl We can see that virtual 
crystal approximation quantitatively reproduces the lattice constants and the energy differ- 
ences between various magnetic orderings when compared to the calculation with "real" La 
and Sr atoms. The accuracy of magnetic energy differences is consistent with earlier work 



23]. 



Appendix B: The effects of SrTiC>3 substrate 

Inclusion of SrTiC>3 substrate in the simulation of PbTi03/Lai_ x Sr x Mn03 heterostruc- 
tures is computationally intensive. In this section, we test the effects of SrTiC>3 substrate 
on structrual and magnetic properties on the PbTiOs/Lai-^Sr^MnOs interface. We com- 
pare two calculations: one with three unit cells of SrTiC>3 and the other without SrTiCV 
Lai_j.Srj.MnO3 is four unit cells thick with nominal doping x = 0.2. PbTiC>3 thin film is 
polarized so that both calculations are in the accumulation state. In terms of structural 
properties, we focus on two important quantities. One is the c/a ratio of each oxygen oc- 
tahedron that encloses Mn atoms. The other is the rumpling 5 of each Mn02 layer. The 
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TABLE VII: Comparison between the virtual crystal approximation and supercell calculations. 
The nominal doping x is 0.5 in the virtual crystal approximation. A c(2 x 2) x 2 supercell is 
employed with La and Sr atoms forming a checker-board pattern (every nearest neighbor of Sr is 
La and vice versa). The lattice constants reported are those for A- type antiferromagnetic ordering. 
AE is the energy difference between ferromagnetic ordering and A-type antiferromagnetic ordering 
per Mn atom, defined by Eq. (pQ). A range of Hubbard U (0 < U < 2 eV) are tested. 





Virtual crystal approximation 


c(2 x 2) x 2 supercell 


U (eV) 


a (A) 


6(A) 


c(A) 


AE (meV) 


a (A) 


6(A) 


c(A) 


AE (meV) 




F 


A 


F 


A 


F 


A 




F 


A 


F 


A 


F 


A 







5.345 


5.363 


7.561 


7.385 


5.386 


5.408 


-34 


5.327 


5.366 


7.482 


7.372 


5.381 


5.411 


-37 


1 


5.385 


5.376 


7.612 


7.411 


5.430 


5.433 


-10 


5.355 


5.380 


7.513 


7.398 


5.400 


5.439 


-11 


2 


5.401 


5.380 


7.623 


7.422 


5.434 


5.452 


7 


5.365 


5.387 


7.520 


7.411 


5.410 


5.452 


6 



results are shown in Fig. ITUA., Layer 1 is the interface. Layer 4 is the artificial surface when 
SrTiC>3 is absent. The solid symbols are with SrTiC>3 and the open symbols are without 
SrTiOs. We can see that the structural properties with or without SrTiOs substrate quickly 
converge as the interface is approached. At Layer 1, the difference is negligible. In terms of 
magnetic properties, we calculate the <i-orbital magnetic moment of each Mn atom, using 
the Lowdin orbitals approach |16 |. Fig. ITTB shows the comparison between with SrTiC>3 
substrate (solid symbols) and without SrTi0 3 substrate (open symbols). Similar to struc- 
tural properties, the effects of SrTiOs substrate on magnetic properties are generally very 
small and diminish at the interface. Therefore we do not include SrTiOs substrate in our 
simulation, not only to reduce the computation burden but also to introduce an artificial 
surface so that we can apply a counting method (see Appendix [U]) to accurately calculate 
hole distribution in Lai_ x Sr x MnC>3. 
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B 





FIG. 11: (Color online) The effects of SrTiO"3 substrate on structural and magnetic properties at 
the PbTi03/Lai_ x Sr x Mn03 interface. Layerl is the interfacial layer of Mn02- The solid symbols 
correspond to the data with SrTi03 substrate. The open symbols correspond to the data without 
SrTi03 substrate. A) rumplings of each Mn02 layer and c/a ratio of each oxygen octahedron that 
encloses Mn atoms. B) ci-orbital magnetic moment of each Mn atom, calculated by using Lowdin 
orbit als. 

Appendix C: The method for counting holes 

A widely used approach for calculating the number of holes and the magnetization of Mn 
atoms is to use Lowdin orbitals 16| . However, a more direct method is to use the electron 



density itself. The difficulty lies in that the boundary between each manganite layer is not 
well-defined in the thin film of Lai-^Sr^MnOs. We develop a method to self-consistently set 
the boundary between each manganite layer, provided that the manganite is half-metallic. 

For a half-metallic manganite, there are no states at the Fermi level in the minority spin 
channel, so that there must be a definite integer number of electrons N c filled in the minority 
spin channel. N c depends on the details of pseudopotentials. For our pseudo atoms (see 
Tabled]), for doping level x, we have La 3+ : 5s 2 5p 6 5<i 6s , Sr 2+ : 4s 2 4p 6 5s°, 2 ~: 2s 2 2p 6 and 
due to charge conservation, Mn ion is nominally +(3 + x) with an electron configuration 
3s 2 3p 6 3d 4 ~ x 4:S . For one unit cell of Lai-^Sr^MnOs, since spin polarization only comes from 
the electrons on Mn <i-orbitals, we can sum all the other electrons that are formally spin 
unpolarized: 8x (1 — x) + 8xx + 8x3 + 8 = 40. The four terms are from La 3+ , Sr 2+ , three 
2 ~ and the Mn ion without ci-electrons. Hence N c = 40/2 = 20, which is independent of 
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FIG. 12: (Color online) Illustration of how to count the charge in each layer of Lai-^Sr^MnOa. The 
spin polarized part is Lai-^Sr^MnOs. The green dashed lines highlight the computed boundaries 
of each layer. 

hole doping. Now we start from the vacuum (see Fig. [T2|) where there is no charge. We 
integrate the minority spin channel moving into the film until the integral is equal to 20. 
Then this position determines the boundary of the first layer. Next we restart the integral 
from this boundary until it reaches 20 again. This determines the boundary of the second 
layer. Repeating the procedure yields the boundaries of each manganite layer. Once the 
boundaries are determined, we integrate the charge density of both majority and minority 
spins in each layer and thus layer-resolved holes and magntization follow straightforwardly. 

We comment that in DFT calculations, as long as the Hubbard U is larger than a critical 
value U c , the manganites become half-metallic in the ferromagnetic phase. For SrTiCv 
strained Pnma Lai-^Sr^MnOa, we find U c ~ 1 eV. Therefore, for the useful and reasonable 
range of U, our method is valid. 
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FIG. 13: (Color online) The complete magnetic phase diagram of Lai-^Sr^MnOs (4 formula, 
20 atoms). The calculation is based on LSDA+C7 with U = 1 eV. AE is the energy difference 
between ferromagnetism and various types of anti-ferromagnetism. The labels F, A, C and G 
refer to ferromagnetism, A-type, C-type and G-type antiferromagnetism, respectively. Each label 
highlights the ground state magnetic structure of the given hole doping region. 

Appendix D: The phase diagram of manganites from LSDA+C7 

In this section, we provide the complete magnetic phase diagram of Lai-^Sr^MnOs 
in Fig. H3] based on LSDA+[7 with U = 1 eV. In the high doping region (x > 0.4), 



LSDA 

states 



U does reproduce the experimentally observed sequence of different magnetic ground 



23|: FM -> A - AFM -> C - AFM ->■ G - AFM as the hole doping x increases. Our 



result is consistent with previous calculations 



12|. 



Appendix E: Test of band alignment and possible artificial charge spillage 

Due to the well known underestimation of band gaps in DFT calculations, band align- 
ment errors and possible artificial charge spillage into the conduction bands of on material 
at an interface may occur. These errors can lead to unrealistic ground states when sim- 



ulating the interface between ferroelectrics and metals 
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37| . We check our calculations of 
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FIG. 14: (Color online) Atomic projected density of states (PDOS) of PbTiO"3 in the 
PbTiOa/Lai-zSr^MOa/Pt heterostructures. a) The accumulation state, b) The depletion state. 
The red curves are Ti-d projected states and the blue curves are O-p projected states. Layerl 
refers to the interface between PbTiOs and Lai_ x Sr x Mn03 (the interface we are interested in) and 
Layer4 is the interface between PbTiOs and Pt. The green solid line is the Fermi level. 

PbTiOs/Lai-zSr^MOa/Pt interface and find that the Fermi level is in the band gap of 
PbTiOs. A typical projected density of states (PDOS) of both the accumulation and de- 
pletion states is illustrated in Fig. [131 where Layerl refers to the interface between PbTiOs 
and Lai-xSr^MnOs (the interface we are interested in) and Layer4 is the interface between 
PbTiOs and Pt. We can see that the interior of PbTiOs remains insulating. We need to 
point out that both terminations of PbTiOs are PbO layers in our calculations, instead of 
the pathological Ti02 termination which leads to a metallic ferroelectric ground state in 
other similar systems 371 ] . 
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